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ABSTRACT 

We present calculations of magnetic field evolution by the Hall effect and Ohmic decay in the crust of neutron 
stars (NSs). In accreting NSs, Ohmic decay is always the dominant effect due to the large resistivity. In isolated 
NSs with relatively pure crusts, the Hall effect dominates Ohmic decay after a time fswitch — 10"* yr B'^^, where 
Bi2 is the magnetic field strength in units of 10'^ G. We compute the evolution of an initial field distribution by 
Ohmic decay, and give approximate analytic formulas for both the surface and interior fields as a function of 
time. Due to the strong dependence of fswitch on B12, early Ohmic decay can alter the currents down to the base 
of the crust for B ^ 10' 'G, neutron drip forZ? 10'^ G, and near the top of the crust for B > 10'^ G, respectively. 
We then discuss magnetic field evolution by the Hall effect. Several examples are given to illustrate how an 
initial field configuration evolves. Hall wave eigenfunctions are computed, including the effect of the large 
density change across the crust. We estimate the response of the crust to the magnetic stresses induced by 
Hall waves, and give a detailed discussion of the boundary conditions at the solid-liquid interface. Finally, we 
discuss the implications for the Hall cascade proposed by Goldreich & Reisenegger. 
Subject headings: stars:magnetic fields — stars:neutron 



1. INTRODUCTION 

Our picture of neutron star magnetic field evolution comes 
from studies of radio pulsars, accreting neutron stars, and 
most recently, magnetars. The known radio pulsars fall into 
two general classes: young (< 10^ yrs) pulsars with spin pe- 
riods ^ 1 s and ^ 10'^ G magnetic fields, and the millisecond 
radio pulsars, which have inferred magnetic fields as low as 
10** G. Whereas most radio pulsars are isolated objects, the 
millisecond pulsars are predominantly in binaries and have 
undergone a period of accretion. This has led to the sugges- 
tion that millisecond pulsars form by accretion, which both 
spins up the neutron start to short rotation periods, and also 
acts to reduce the magnetic field (see Bhattacharya 1995 for a 
review). 

This "recycling" scenario is consistent with observations 
of accreting neutron stars. Young accreting neutron stars in 
high mass X-ray binaries (HMXBs) are observed as X-ray 
pulsars with inferred 10'^ G magnetic fields, confirmed by 
detection of cyclotron lines in their X-ray spectra (Triimper 
et al. 1978). Old accreting neutron stars in low mass X-ray 
binaries (LMXBs) are generally not X-ray pulsars, implying 
no magnetic disruption of the accretion disk in those systems, 
and magnetic fields < 10^-10'" G (see Gumming, Zweibel, 
& Bildsten 2001 for a recent discussion). Recently, 5 mil- 
lisecond accreting X-ray pulsars have been discovered in tran- 
sient LMXBs, with spin periods between 180 and 410 Hz, 
and inferred magnetic fields 10^ — lO' G, exactly in the range 
expected for millisecond pulsar progenitors (Chakrabarty et 
al. 2003). 

An additional possibility is that neutron star magnetic fields 
decay with age. Field evolution in isolated radio pulsars has 
been studied using their space and velocity distributions, as 



well as the statistical distribution in the P-P diagram. This 
probes field evolution on the timescale of the ages of radio 
pulsars (< 10^ yrs). At first, it was argued that exponential 
field decay on a few million year timescale was required to un- 
derstand the P-P distribution (Ostriker & Gunn 1969; Gunn 
& Ostriker 1970; Narayan & Ostriker 1990), but the work of 
Bhattacharya et al. (1992) and others in the 90's argued that 
in fact no field decay occurs over the lifetime of a radio pul- 
sar. However, this issue remains somewhat open: recently, 
several authors have again raised the suggestion that neutron 
star magnetic fields evolve on a timescale something like 10 
million years (Cordes & Chernoff 1998; Tauris & Manchester 
1998; Tauris & Konai- 2001; Gonthier et al. 2002). 

The discovery of magnetars has complicated this picture, 
but also provided a new opportunity to study magnetic field 
evolution on relatively short timescales. Observed as the 
young < lO'* yr soft gamma repeaters (SGRs) and anoma- 
lous X-ray pulsars (AXPs), these sources are believed to be 
directly powered by decay of ultrastrong ^ 10''*-10'^ G mag- 
netic fields. Long term release of magnetic energy may act to 
delay the cooling of the neutron star (Thompson & Duncan 
1996; Heyl & Kulkai-ni 1998; Colpi, Geppert, & Page 2000; 
Arras, Gumming, and Thompson 2004), while rapid release of 
magnetic stresses that build in the neutron star crust (Duncan 
& Thompson 1992; Thompson & Duncan 1995, 1996) may 
power the frequent gamma-ray bursts observed from SGRs, 
and now also AXP IE 1048.1-5937 (Gavril, Kaspi, & Woods 
2002). 

These ideas and observations have motivated a large range 
of theoretical work on neutron star magnetic field evolution 
(e.g. see Bhattacharya & Srinivasan 1995 for a review). In 
this paper, we discuss the evolution of currents in the neutron 
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star crust, which evolve due to Ohmic decay, and by the Hall 
effect. The importance of currents in the crust is that they de- 
cay much more quickly than currents in the core. They there- 
fore play an important role in models of accretion-induced 
field evolution: either the currents are assumed to flow only 
within in the crust, for example produced after the birth of 
the neutron star by thermomagnetic effects (Blandford, Ap- 
plegate, & Hernquist 1983), or magnetic flux may be expelled 
from the core into the crust due to vortex-fluxoid interactions 
(e.g. Srinivasan et al. 1990; Ruderman, Zhu, & Chen 1998). 
At higher field strengths. Hall effects most likely are respon- 
sible for building up stresses in magnetar crusts which then 
are released during SGR and AXP bursts. 

Previous investigations into Ohmic decay have focused 
mainly on the decay of the surface field, and involved cal- 
culating Ohmic decay eigenmodes for the crust (Sang & 
Chanmugam 1987), self-similar solutions near the surface of 
the star (Urpin, Chanmugam, & Sang 1994), or direct time- 
integrations of the Ohmic diffusion equation for an initial dis- 
tribution of currents (Urpin & Muslimov 1992; Urpin & Ko- 
nenkov 1997; Page, Geppert, & Zannias 2000). For accreting 
stars, the evolution due to Ohmic decay for fields supported 
by crustal currents has been calculated in great detail (e.g., 
Geppert & Urpin 1994; Urpin, Geppert, & Konenkov 1998; 
Konar & Bhattacharya 1997). We revisit Ohmic decay in this 
paper, calculating both Ohmic decay modes and self-similar 
solutions, and show that they are related in a simple way. As 
well as constraining the "initial condition" for the Hall effect, 
this allows a general way to understand the distribution of cur- 
rents in the crust at a given age, as well as the decline of the 
surface field. 

The evolution by the Hall effect is much less well- 
understood, being non-linear and therefore much more com- 
plex. Although the Hall effect itself is non-dissipative, Gol- 
dreich & Reisenegger (1992, hereafter GR) proposed that it 
leads to dissipation through a turbulent "Hall cascade", mag- 
netic energy cascading from large to small scales where it 
dissipates by Ohmic decay in the Hall time associated with 
the largest lengthscale. Recent studies by Biskamp and col- 
laborators (Biskamp et al. 1999) of Whistler turbulence in- 
deed show a clear cascade of energy over more than an or- 
der of magnitude in lengthscale. Other authors have also 
demonstrated transfer of magnetic energy between different 
scales (Muslimov 1994; Naito & Kojima 1994; Shalybkov & 
Urpin 1997; Urpin & Shalybkov 1999; Hollerbach & Rudiger 
2002). In a different approach, Vainshtein, Chitre, & Olinto 
(2000) included the effect of the density gradient in a simple 
model. Most recently, Rheinhardt & Geppert (2000), Geppert 
& Rheinhardt (2002), and Rheinhardt, Konenkov, & Geppert 
(2003) discuss a "Hall drift instability" which can lead to non- 
local transfer of energy to small scales. 

In this paper, we study the implications of realistic crust 
models for evolution of crustal currents. In §2, we calculate 
the timescales associated with Ohmic decay and the Hall ef- 
fect for realistic crust models. We discuss the relative impor- 
tance of Ohmic decay and the Hall effect at different stages of 
a neutron star's life. In §3, we consider Ohmic decay in detail. 
We show that previous studies of Ohmic decay can be under- 
stood in a simple way. In §4, we discuss the evolution due to 
the Hall effect. We start with a simple discussion of the initial 
evolution of a dipole field which gives a simple physical pic- 
ture. We then calculate the properties of linear Hall waves in 
a realistic crust model. We clarify the physical nature of the 
"Hall drift instability" and comment on its relevance for the 



overall field evolution. We summarize our main results in §5. 
In the Appendix, we discuss the behavior of Hall waves at the 
boundaries of the crust. 

2. TIMESCALES FOR THE HALL EFFECT AND OHMIC DECAY 

In this section, we estimate the timescales of evolution by 
the Hall effect and Ohmic decay in both isolated and accreting 
neutron stars, and discuss the regimes where the Hall effect 
is important. Both the Hall and Ohmic timescales are most 
sensitive to the lengthscale L over which the magnetic field 
and crust properties vary. In this section we make a simple 
estimate by setting L equal to the local pressure scale height 
H. This allows a comparison between the two processes, and 
sets the scene for the more detailed discussion in ^and|3 

2.1. Calculation of the Hall timescale 

The evolution of the magnetic field is given by Faraday's 
law, 

dB _ 

— =-cVxE, (1) 
ot 



together with Ohm's law 



^ 1 ^ J JxB 

E = — V X B+ — + , 

c a n^ec 



(2) 



where cr is the electrical conductivity, and the current density 
is J = (c/47r)( V X B) (see GR for a derivation). The first term 
in equation (|3 describes advection of the field by the fluid, the 
second term describes Ohmic dissipation, and the third term is 
the Hall electric field'. Neglecting the advection and Ohmic 
terms, the magnetic field evolution is given by 



dB „ fJxB 
— =-Vx 

at \ UeC 



with a timescale 



fHall ■ 



J 



cB 



(3) 



(4) 



where L is the typical lengthscale over which B, J, and n^. vary. 

The basic physics of equation Q is that the magnetic field 
is frozen to the electrons. In ideal MHD, the electron and ion 
velocities are the same - both follow the E x B drift. In Hall 
MHD, the ions follow separate dynamics. In a plasma or liq- 
uid, this only happens on scales of order the ion skindepth 
(5, = va/uJci (or, in a high [3 plasma, scales of order the ion 
gyroradius Vj/ujci)- In solids, the ions are immobilized by 
lattice forces and Hall effects are important on much larger 
scales. This is the case in the neutron star crust, which elas- 
tically adjusts as the field evolves, balancing the JxB forces 
with shear stresses. We discuss the issue of force balance in 
more detail in §4, where we consider the behavior of wave- 
like solutions at the boundary between the solid crust and fluid 
core/ocean. This is illustrated by Figure[2which shows a vari- 
ant of the usual laboratory demonstration of the Hall effect. A 
current-carrying conductor is placed in an external magnetic 
field, with the current perpendicular to the field. The conduct- 
ing electrons, which move with velocity Ve = -J /neC, are de- 
flected sideways by the magnetic field, and a Hall electric field 
Eh = Ve X B/c develops to balance the magnetic force. The 

' We do not discuss the battery or thermoelectric terms in this paper 
(e.g. Blandford et al. 1983) since they are suppressed by a factor {kgT /Ef j,)^ 
for a degenerate gas, and are important in the crust only for very young neu- 
tron stars. 
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Fig. 1. — Simple example of the Hall effect: (a) the usual laboratory 
set up in which the Hall electric field balances the magnetic force on the 
conducting electrons; (b) if the conductivity varies with height, the elec- 
tron velocity depends on height, shearing the magnetic field on a timescale 
'Hall = L/ve = neeL/J. 



variation on the usual experiment is to allow the conductivity 
to depend on height. In this case, J and therefore v^. and Eh 
depend on height, giving an emf S = j En-dl which causes a 
magnetic field component to grow along the current direction. 
The result is that the magnetic field lines are sheared by the 
electron flow. In the star, this shear arises because of gradients 
in B, J, and Ug. The Hall timescale given in equation @ is 
simply the shearing time of the electrons, feaii = I^/ve- 

We now calculate fnaii for densities p from « 10^ " o„i-3 
through neutron drip at p r 
the crust/core interface at p 



gem 

4 X 10" g cm"^, down to 

< 



10"* gcm-l 



For 



6 X 10'^ g cm ^, the relativistic, degenerate electrons dom- 
inate the pressure. The electron Fermi energy and 

1/3 1/3 

pressure are Ef ^, = 51 MeV Ye , and Pe = 1.2 x 



10 erg cm 



where pn = p/lO'^ gem and Ye 



is the number fraction of electrons. The top of the crust is set 
by r = Z^e^/akBT = 175 (Potekhin & Chabrier 2000), where 
a = (3/4™,)'/'' is the interion spacing, giving a density 



X 10^ 



gem 



26 



(5) 



The top of the crust varies strongly with temperature, from 
p ~ 10^ g cm"-' in accreting stars, to 10^ g cm"-' in isolated, 
cooling neutron stars. 

For /O^6xl0'~g cm"-', the degenerate, non-relativistic 
neutrons dominate the pressure. The neutron Fermi energy 

and pressure are Ef^„ = 15 MeV pl!^'Y^^^(f /0.5), and P„ = 
6.0 X 10^2 erg cm"^ p%^Y,^^\f/0.5). Here, / is a factor which 
accounts for the interactions between neutrons. We use the fit 
of Mackie & Baym (1977) for Ef_„ (following Brown 2000). 

The electron and neutron fractions (Ye and Y„) are given 
by the composition at each depth. An accreting neutron star 
can replace its crust with accreted material, in which case the 



composition is quite different from an isolated object. Fig- 
ure 13 shows the pressure. Ye, Y„ and nuclear charge Z as a 
function of depth in both cases. For an isolated crust, we use 
the cold-catalysed matter results of Haensel & Pichon (1994) 
(p < 4 X 10" g cm"-^), and Douchin & Haensel (2001) (p > 
4 X 10" g cm"-'). For the accreted case, we use the results of 
Haensel & Zdunik (1990, hereafter HZ), who calculated the 
composition changes due to electron captures and pycnonu- 
clear reactions as fluid elements traverse the crust during ac- 
cretion. HZ stop their calculation at p= 1.25 x 10'-' g cm"-'; 
at higher densities, we fix the composition at their final value, 
giving Ye = 0.05 for p > 1.25 x 10'^ g cm"^ in our accreted 
crust. Our equation of state agrees well with Negele & Vau- 
therin(1973). 

We show the Hall timescale as a function of density for an 
isolated and accreted crust in Figure|3 We choose the length- 
scale L to be the local pressure scale height, H = P/ pg, where 
g is the local gravity, assumed constant. For densities between 
10'~ and 10'-' g cm"-', the Hall time is approximately constant 
with increasing density because as neutron drip occurs, the 
electron fraction drops dramatically, cancelling the increasing 
density and scale height. At lower densities, the electrons set 

tiie pressure, giving H = 77.6 m p[{\Ye / 0. 25)'^ /\2A5/g 14), 
and 



_ 5.7 X lO^yrs 5/3 / Y, 

^Hall, outer — ~ P\2 



B 



12 



0.25 



11/3 



2.45 



(6) 



At densities greater than neutron drip, the neutrons dominate 
the pressure, giving 

H: 



: 245 m p%\f/0.5)Y^^\2A5/gi4), and 



^Hall, inner ~ 



1.2 X 10^ yrs 7/3^10/3 
on 



0.05 



.2.45 



(7) 

We have inserted typical values for Ye at each density. In ^ 
we argue that the timescale for overall evolution of the field 
is roughly the Hall time at the base of the crust, as given by 
equation 0. The choice L = H introduces the 1/g depen- 
dence of fHaii on the local gravity in the crust, and therefore 
some sensitivity to the equation of state of the core. 

2.2. Calculation of the Ohmic decay time 

For comparison, we now calculate the Ohmic time across 
a scale height fohm = ^ttcfH^ /c^. The conductivity is cr = 

HeC- /m^,v, or 



cr=1.52x 10^^ s"' 



10'7 s"i 



2/3 2/3 



Pn 



(8) 



where v is the electron scattering frequency, and = 
Ef/c^ 3> nie is the electron effective mass, giving fohm = 
4x10^ yrs fxis^^l- 

The electrical conductivity is set by electron scattering from 
phonons and impurities. The phonon scattering frequency is 
Uep = (l3akBTAep/n)[l+(e/3.5Tf]-^^^exp(-Tu/T) (Baiko 
& Yakovlev 1995, 1996), where Aep is the order unity 
Coulomb logarithm, and is the Debye temperature. 



0.45/1 MTrZ^e^n,. 



1/2 



= 3.5 X 10''KF,p[f , (9) 



with Aceii the number of baryons per unit cell. Except for 
very young neutron stars, T ^ Qo- The exponential sup- 
pression of the scattering frequency for T < Tu = %.! x 
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Fig. 2. — Pressure, neutron fraction Y„, nuclear charge Z, and electron frac- 
tion Ff. as a function of density for an isolated crust (solid line) and accreted 
crust (dotted line). 



10^ Kp[f (};/0.05)(Z/30)'/^ occurs because Umklapp pro- 
cesses freeze out below this temperature. For T > Tu, we find 

f 1.2 X 10^4 ,-i ipy'/Ti)(Yj0.25)'/' 
' \l.8xl02^s-i(p;/Vr,2)(V0.05)^/3- 

The impurity scattering frequency is (Itoh & Kohyama 
1993), iy,Q = 1.75 X lO"' s"' xQA,q/Z, where x = pr/m^c = 

100 pI^^I^^ , is the order unity Coulomb logarithm, and 
Q is the impurity factor, giving 



^ ^ f 1.6 X 1025 {p\l^'/Q)(Ye/Q.25fl\Z/3Q) 
^ I 4.4 X 1025 {p\l^ /Q)(Ye/0.05yl\Z/3Q) 



(11) 



For clarity, we set K^p = A^q = 1 in equations ( I10> and (II l\ : in 
our numerical calculations, we calculate these factors accord- 
ing to Itoh & Kohyama (1993) and Baiko & Yakovlev (1996), 
and include the exponential suppression of phonon scattering 
for r < 71/ as in Gnedin, Yakovlev, & Potekhin (2001). Our 
conductivities agree well with those of Potekhin et al. (1999). 



Impurity scattering dominates either at low temperatures 
T < Tif (when phonon scattering is suppressed), or for high 
impurity levels Q > Qait, where 

2.4 (Ti/p%')(Z/30)(0.05 
\l3(Ti/pli')(Z/30){0.25/Y,r'' 

is given by the ratio of equations jlO> and il 1> . For an equi- 
librium crust in an isolated neutron star. Flowers & Ruderman 
(1977) estimated Q w 10"^, in which case the conductivity is 
set by phonon scattering above 10^"** K, depending on posi- 
tion in the crust, and by impurity scattering at lower tempera- 
tures. However, recent calculations by Jones (2001) suggest a 
much larger value Q ^ 10, in which case impurity scattering 
may set the conductivity even for T > 10^ K. 

In accreting neutron stars, the crust temperature is high 
(> 10** K) for rapid accretors (> 10"" M© yr"'), in which 
case phonons dominate for Q ^ 0{l). However, Q is set by 
the composition of the nuclear burning which occurs at low 
densities, and is likely to be > (9(1) (for example, Schatz et 
al. 1999 found Q « 100), in which case impurities dominate. 

When impurities dominate, the Ohmic time is 

0.05 j V 0.8 



ro^ = 5.7Myrs^(^- 



0.5 / V2.45/ ' 



and when phonons dominate, it is 



15/6 



fOhn 



Y \ 5/3 / y \ 10/3 



:2.2Myrs^ ^ [ 

^ Ti I 0.05/ I 0.8 



Xrflll)-^ (14) 
0.5/ V2.45/ ^ ' 



(15) 



The ohmic time fohm = 4-TTaL?/c^ is shown in Figure |3] for 
equilibrium crusts with g = 10"^ and T = 10^ K and 10** K, 
and for an accreted crust with T = 5 x 10** K and Q=l. 

2.3. When does the Hall Effect dominate Ohmic decay? 

The ratio of Ohmic time fohm to the Hall time t^ is equiva- 
lent to the product of electron cyclotron frequency and elec- 
tron collision time, fir = eBr/mi^c. Large ilr means that 
fHaii <C ?ohm, and the Hall effect dominates. If the conductivity 
is set by impurities {Q > Qcnt 01 T <Tu) 

. _ / 7 (Bn/fip'f )(Z/30)(};/0.25)-2/3 
. 1 (Bn/epJ )(Z/30)(y,/0.05)-2/3 ' 

and if the conductivity is set by phonons (Q < Q^n and T > 

Tu) 

, _ / 0.6 (Bnp/2/Ti)(Ye/0.25f/^exp(Tu/T) 
0.4 (Bnp\i'/Ti)(Ye/0.05)y'exp(Tu/T) ' 

These different regimes are summarized in Figure |3 We 
show the fir = 1 contour in the B-T plane for conditions near 
the base of the crust (p = 10'"* g cm"-'). Note, however, that 
rtr is reasonably insensitive to p, so that all depths change 
from Ohmic decay-dominated to Hall-dominated at roughly 
the same time. 

Figure 0] shows very clearly the difference in evolution be- 
tween isolated and accreting neutron stars. An isolated neu- 
tron star cools to T < 7[/ in w 1 Myr, after which impurities 
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Fig. 3. — Hall timescale fnaii (left panel) and ohmic time ?(,hm (right panel) as a function of density for two examples of an equilibrium crust (solid lines) 
and accreted crust (dashed line). The equilibrium crusts have Q — 10"^ and T — 10* K and 10^ K; the accreted crust has T = 5 X 10^ K and 2=1. We assume 
the crust is isothermal, with a constant magnetic field B = 10'^ G ((Hall ^/B), and take the relevant lengthscale to be the scale height, with local gravity 
g = 2.45 X 10'** cm s^^. 
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Fig. 4. — Ohmic decay vs. Hall effect. The solid lines show the B 
field above which Qt > 1 as a function of temperature at a density of 
p = lO'* g cm"', and for an isolated crust. 



set the conductivity giving fir 2± B12/Q at the base. The ex- 
tent to which Hall effects dominate then depends sensitively 
on the value of Q for an isolated crust. The recent calcu- 
lation by Jones (2001) of the formation enthalpies of point 
defects in the crust gives 2 > 1 at densities greater than neu- 
tron drip (see also de Blasio 2000 for densities lower than 
neutron drip). In this case Ohmic decay dominates the evolu- 
tion even at late times. The original estimates of Flowers & 
Ruderman (1977) gave Q in which case the Hall effect 



dominates for T <Tu- An accreting neutron star has T > Tu 
because of heating by accretion (e.g. Brown 2000), and in ad- 
dition may have a large impurity parameter Q> I because of 
the complex mixture of nuclei produced by burning of the ac- 
creted hydrogen and helium (Schatz et al. 1999). This means 
that whether impurities or phonons dominate, the conductiv- 
ity is poor enough that Ht <C 1, implying that the Hall effect 
is not important for an accreting star (even after accretion ter- 
minates, if the crust is impure). 

For an isolated neutron star with a pure crust (2 <C 1), Fig- 
ure 0] allows us to estimate when the Hall effect will dom- 
inate Ohmic decay. For Q = 10"^ at the base of the crust, 
the critical temperature below which the Hall effect is domi- 
nant is r « 6 X 10^,2 X 10*^, and 6 x lO^K for magnetic fields 
B = 10", 10'^ and lO'"* G, respectively. At these high temper- 
atures, T >Tu and the conductivity is set by electron-phonon 
scattering. To estimate the core temperature as a function of 
time, we adopt the cooling law for modified URCA neutrino 
emission with specific heat given by normal (unpaired) neu- 
trons, t ~ 10'' yr T^^ (e.g. Shapiro and Teukolsky 1983). Us- 
ing equation ( I16t . the age at which the Hall effect becomes 
dominant is 

fswitch ^2 X 10^1-^12 • (17) 

Ohmic decay is more important in lower field stars since the 
Hall effect takes over at a later time. For typical pulsar mag- 
netic fields, the Hall effect starts to dominate while the star is 
still quite young. Ohmic decay calculations ignoring the Hall 
effect are only vaUd for t < f switch ■ 

3. FIELD EVOLUTION WHEN Qt < 1 

In this section, we compute the Ohmic decay of currents in 
the crust at early times when Ht <C 1, before the Hall effect 
dominates the field evolution. We first calculate Ohmic decay 
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modes, and then discuss the Green's function solution, which 
shows self-similar behavior at late times. We show that these 
solutions are simply related, giving a general way to under- 
stand the distribution of currents in the crust at a given age, 
as well as the decline of the surface field. This constrains the 
"initial condition" for the later evolution by the Hall effect. 

3.1. Ohmic decay modes 

We now calculate the Ohmic decay modes for currents con- 
fined to the crust. We write the poloidal magnetic field in 
terms of the potential /(r,f), S = Vx Vx [/(r,f)}'"'(6',0)] 
(Chandrasekhar 1961). The field components are 



m+i) 



fYrer + 



Idf fdYl" dY; 



60 + - 



r dr \ 86 " dcf) 
Faraday's law (eq. (Q) with E = J/a gives 



5/ 
dt 



= r]{r,t) 



Ul + \) 



(18) 



(19) 



where rj = /Ana is the magnetic diffusivity. 

When Ohmic decay determines the evolution of the field, 
the conductivity is mainly set by phonons, in which case 
the time-dependence of rj arises only through the temper- 
ature, rj oc T^. Following Urpin et al. (1994), we absorb 
this time-dependence into a new time variable f, given by 
di ex dt. For modified URCA cooling T oc f"'/^, we find 
f = (f/1.5 Myr)-'/^ Myr when 77 is evaluated at the reference 
temperature Tg = 1. We then look for Ohmic decay modes, 
/ oc exp(-f7r), giving 



d^f _Kl+\) 
'dr^~~ 



f = 



rjT 



(20) 



Note that exponential decay in i represents a slower decay 
in real time f, because the conductivity of the cooling crust 
grows with time, reducing the Ohmic dissipation rate. 

Equation (I20> is an eigenvalue problem for t given bound- 
ary conditions for /. At the top of the crust, we match onto 
a vacuum field, giving / oc r"', and normalize the field such 
that Br = 1 there. At the base of the crust, we consider two 
different boundary conditions, either = at the base, im- 
plying the field is confined entirely to the crust, or we assume 
the field connects to an interior vacuum field, an example of a 
field which threads the neutron star core^. 

A WKB analysis clarifies the basic properties of the modes. 
Writing / oc exp{i drkr)exp(-t /t), we find the dispersion 
relation 

*; = _L-fl^.„.,-. ,21, 

rjT 

where k,- is the radial wavenumber and r is the decay time. 
We ignore the horizontal wavenumber term, which is negli- 
gible unless I > R/Hh ^ 20. The WKB solution, including 
amplitude and phase, is 



/ oc ry'/'^exp ( ±/ 



dr 



(rjT) 



1/2 



(22) 



- In fact, if the field penetrates the core, then there is no strict boundary 
condition quantizing the phase in units of roughly nn in the crust. The vertical 
wavelength in the core will be much smaller than in the crust, and then the 
value and slope at the crust core interface will reflect quantizing the phase in 
the core. For eigenmodes with many nodes in the crust, this will not change 
the waves significantly, but for long decay times 3> 10* yr, this boundary 
condition is crucial. 



Table 1. Ohmic Decay Modes (dipole) 



fl 


'Ohm/T„ 


'Ohm/'l Tn 




= 1,2= 10^3 


tnu„, = 7 5 X 10* vrs 


1 


3.60 


3.60 


2 


18.2 


4.56 


3 


44.0 


4.89 


4 


82.0 


5.12 


5 


133.0 


5.33 


7i 


= 3, Q = 10-'' 


tnu,„ = 7 X 10' vrs 


1 


2.32 


2.32 


2 


12.1 


3.03 


3 


29.8 


3.31 


4 


55.8 


3.49 


5 


90.4 


3.61 




:0.oi, g= ur 


^'0hm=4.7x lO^Jyrs 




0.315 


0.315 


2 


2.17 


0.542 


3 


5.67 


0.630 


4 


10.8 


0.675 


5 


17.6 


0.703 




= 1,2=10-', 


, ?ohm = V-5 X 10" yrs 




inner dipole boundary condition 


1 


0.0729 


0.0729 


2 


8.76 


2.19 


3 


28.8 


3.20 


4 


60.5 


3.78 


5 


105 


4.20 



For /, the amplitude and wavelength grow oc 77^/^ toward the 
top of the crust. The magnetic field and current components 
have the same oscillatory phase, but their envelopes scale as 
B,- oc 77'/'*, Z?e,0 oc T?"'/"^, and Jg,^ oc ri~^^*. A WKB turning 
point is reached at the depth Ztpir) at which kMa ^ 1, or 
equivalently, where the mode decay time equals the local dif- 
fusion time. Above this depth, the left and right sides of equa- 
tion ilOi can no longer balance, leading to evanescent behav- 
ior 

We show two examples of Ohmic decay eigenmodes with 
Z = 1 in Figures |5l (B,- = at the base) and |6l (vacuum field at 
the base). The left panel in each case shows the lowest order 
eigenmode. When B, = at the base, Br varies rapidly with 
depth, giving Bg w {R/H)Br » B,- at the base; for a vacuum 
boundary condition at the base, B, varies only slightly across 
the crust, giving B, ^ Bg at the base. For short wavelength 
modes (right panels), the eigenfunctions become insensitive 
to the lower boundary condition; the effect of the boundary 
condition is to introduce an overall phase shift. The WKB 
amplitude scalings and turning point can be seen clearly for 
these modes. 

Table [T] gives Ohmic decay times for numbers of nodes 
n = 1-5, different temperatures, impurity fractions and in- 
ner boundary conditions. The results are scaled in terms of 
the Ohmic time at the base, which is listed in the table. For 
T = 10^ K, the decay times for the lowest order modes are 
Ti = 2. 1 X 10'' yr and ri = 1 .0 x 10** yr for the B,. = and vac- 
uum inner boundary conditions. As pointed out by Pethick & 
Sahrling (1995), the difference between these two cases is that 
the lowest order mode has a decay time w R/H times larger 
when the field penetrates into the core, because the length- 
scale on which B, changes through the crust is the radius in- 
stead of the crust thickness. 

3.2. Green's function solution 
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Fig. 5. — Ohmic decay mode eigenfunctions for an equilibrium crust with T = K and g = 10 ^ . We show the lowest order mode {n = 1) in the left panel 
and the n = 5 mode in the right panel. The inner boundary condition is = 0. 




An alternative method of solving the diffusion equation is to law diffusivity rj cx z"" (where z is the depth). ^. We start by 
find the Green's function. We now calculate the eigenmodes 

and Green's function for a plane parallel crust with power- ^ Although we choose a power law diffusivity mainly because it allows an 

analytic solution, this is also a good approximation since the conductivity is 
a power law in both position and time over a large range of scales. 
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Fig. 7. — Green's function for Ohmic decay of poloidal magnetic field. 
Tlie conductivity profile is a = 7/2. Each line is labelled with the value of 



defining the similarity variable 

Jo 



dz 



(23) 



Useful quantities are fdiff(z) = £, t oc z , the characteristic 
diffusion time to depth z, and Zdis(t) oc f jj^g ^^pi]^ fgj- 
which the diffusion time is the current age, i = fjiff, given by 

C(Zdiff,f)= 1- 

Setting ?7 oc z we look for solutions oc /(^) exp(-f /t), 
where S, = S,(z, r), giving 

./(o=r-/-.(a (24) 

where v = (a + 2)"', and J-^ is a Bessel function of the first 
kind. The solution with order -ly becomes constant near the 
surface, therefore matching onto a vacuum field at low den- 
sities. The decay time r for each mode is determined by en- 
forcing the boundary condition at the base. These eigenmodes 
have similar properties to those determined in §3.1. The turn- 
ing point for each mode occurs at ^ '--^ 1, or depth Zdis{T). 
Since the amplitude of J^iO °^ r'/^ for e > 1, we find 



fccz 



-a/4 



1/4 



at large depths, as expected from the WKB 



scalings (eq. 1221 '). 

The Green's function is obtained by writing an initial delta 
function as a sum over the eigenmodes given by equation 
J24> . and using the Fourier-Bessel inversion formula (Jackson 
1975, pg. 100) and a standard integral over Bessel functions 
(Watson 1966, pg. 395). We find the Green's function solution 
(see also Eichler & Cheng 1989) 

M,t) = j d\n£,iG{W)Mzi) (25) 



where 



(26) 



£, = £,iz,i), £,/ = £,{z/,i), Ii, is a modified Bessel function, and 
/o(z) is the initial field profile. For ^ -C 1, G becomes inde- 
pendent of ^, 



,x 2(1-1/) 



exp 1 -— 



(27) 



Figure0shows the Green's function for a = 7/2 and = 2/ 1 1 , 
appropriate for electron-phonon conductivity above neutron 
drip. 

The evolution in time can be understood by taking various 
limits of equation (I25> . Consider the case where the currents 
extend down to a finite depth Zinit initially, /o(z > Zinit) = 0. 
For early times (or sufficiently large depths), the field is un- 
changed: G« (5(^-0 fori^ ^ 1, giving 



/(z,f)~/o(z) 



C>l,Z>Zdiff(f) 



(28) 



For late times (^ ^ 1), equation Ml\ applies, and the solution 
is independent of depth, with a time-dependence determined 
by the integral over When Zdiff(f) < Zinit. the peak of the 
Green's function at ^' w 1 gives 



/(z,f)-/o(Zdiff(f)) 



C< 1, zdiff(d < 



(29) 



so that the solution is determined by the initial value of the 
field at the current diffusion depth Zdiff(r)- If /o varies slowly 
with depth, then / will vary slowly with time during this 
phase, particularly since Zdiff is a weak power of i. When 
Zdisit) > Zinit, however, the integrand drops to zero at a value 
of ^' < 1, so that equation ( I27> gives G oc (^')^*'"'- Writing 
the solution as an integral over z', the time-dependence factors 
out as 

fiz, i) oc oc r("+i)/("+2) e « 1 , zmit) > Zinit. (30) 

Therefore at late times, the field decays as a power law in time, 
independent of the initial field profile (Sang & Chanmugam 
1987). Since i oc we find that the surface field decays as 



m 

Bq 



fdiff(Zinit) 



(31) 



with 6 = 0.55 {a = 7/2) for Zinit above neutron drip, and 
S = 0.49 (a = 7/4) for Zinit below neutron drip. For complete- 
ness, we also note that when impurity scattering dominates 
the conductivity, f oc f and a ~ 4/3, giving S = 0.7. 

The connection between the eigenmode and Green's func- 
tion solutions is now easily seen. The Green's function solu- 
tion implies that / is constant for z < Zdisii), and unchanged 
for z > Zdiff(d- In terms of Ohmic modes, at a given time f, all 
modes with decay times t which propagate closer to the 
surface than Zdiff(f); have exponentially decayed away. There- 
fore for z < Zdiff(f) the solution is dominated by the mode with 
T = i which has turning point at Zdiff(T). The power law decay 
of the surface field at late times arises because the decay time 
of the mode which dominates the surface field increases with 
time (Sang & Chanmugam 1987). For z > Zdiff(r), the field 
profile is a superposition of all WKB modes with decay times 
T > i. 

3.3. Toroidal fields 

Toroidal fields behave in a qualitatively similar way to 
poloidal fields, but with some important quantitative differ- 
ences. Defining the toroidal field as i? = V X (§(r)F/"er) 
(Chandrasekhar 1961), we find the diffusion equation is 



dg 
dt 



d_ 
dr 



dg 
^ dr 



-V- 



1(1+1) 



with WKB solution 

g oc 77"'''"* exp I ±/ 



dz 
(r7r)i/2 



(32) 



(33) 
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Fig. 8. — Green's function for Ohmic decay of toroidal magnetic field. 
The conductivity profile is a = 7/2. Each line is labelled with the value of 

giving Bg^rj, oc 77"'/^, Jg_^ oc r]~^^^, and Jr cx 7/"'/"^. These scal- 
ings are similar to the poloidal case. For diffusivity 77 oc z"", 
the power-law solution has the same form as in equation 
i24i . but with v = (a + l)/{a + 2), and the Bessel function 
must have order +1^ rather than -1/ to ensure zero field at the 
boundary. The toroidal Green's function is plotted in Fig- 
ure |S] The major difference from the poloidal case is for 
late times, for which the field decays oc f-(2"+3)/(Q+2) ^ f-\.2 
(for a = 7/2), much faster than the poloidal case (compare 
eq. II30I ). In addition, the field decreases rapidly towards the 
surface, g oc ^^"^ oc z""""' , implying the interior toroidal fields 
are much larger than those near the surface. 

3.4. Summary of Ohmic decay 

The results of this section give a very simple way to under- 
stand the distribution of currents and the decay of the surface 
dipole field given the age of the star. For a general conductiv- 
ity profile rjiz), the rescaled time corresponding to the upper 
turning point is 

where z = is at the top of the crust. We plot the physical time 

-3/2 

fdiff OC fjjfj for a detailed crust model with electron-phonon 
scattering in Figure|9] This plot can be read as follows: since 
we have included the fact that the temperature is decreasing in 
time, set the ordinate fjjiff equal to the age f , and then read off 
the density at which z = Zdiff- At lower densities, the currents 
have decayed according to equation (13 1> . and we expect / ap- 
proximately constant, dominated by the Ohmic decay mode 
with T K, t; at higher densities, the initial field has not yet 
evolved significantly, and is dominated by modes with r > f . 
We emphasise again that we assume phonons dominate the 
conductivity in Figure |9l This is a good assumption even if 
the crust is impure 2 > 1, except at the very base of the crust 
(see eq. 1121 ). After « 1 Myr, the crust conductivity becomes 
impurity scattering dominated, and the field decay slows dra- 
matically if g < 1. 

Our results compare well with the time-dependent calcula- 
tions of Urpin and collaborators (Urpin & Muslimov 1992; 




10'° 10" 10'= 1013 

p(g cm-'') 



Fig. 9. — Ohmic time corresponding to the upper turning point as a func- 
tion of distance from the top of the crust. We assume phonon conductivity 
and modified URCA cooling {T oc r"'/*) (for f ^ 1 Myr, impurity scattering 
dominates the conductivity, in which case the Ohmic evolution slows consid- 
erably for Q 1). (i) For a given age t, this plot shows the density at which 
cun'ents are just beginning to evolve due to Ohmic diffusion. At lower densi- 
ties, the field is roughly constant with depth, dominated by the Ohmic decay 
mode with t ^1: aX higher densities, the field has not yet evolved signifi- 
cantly, and is dominated by Ohmic decay modes with t > t. (ii) For currents 
initially distiibuted up to a density pinit, this plot gives the time at which the 
surface field evolution changes from roughly constant with time to a power 
law oc r"-^'. 

Urpin & Konenkov 1997; Page et al. 2000). Their initial 
profile /o(z) varies slowly with depth in the crust, giving 
very little decay of the field for t < fdiff(Zinit) (eq- 1291 ). For 
t > rdiff(zinit), the dccay of the observable dipole field agrees 
well with equation i31\ . We do not reproduce the much longer 
(^ 10^ year) decay timescales found by Chamugam & Sang 
(1987); the reason for this discrepancy is not clear 

Att = fswitch (eq- 1171 ') when the Hall effect begins to dom- 
inate Ohmic decay, we find that the currents must reside at 
a density > 10'" g cm"^ (towards the top of the crust) for 
B = 10'^ G, > 10'2 g cm-3 (near neutron drip) for B = 10'^ G, 
and > lO'** g cm~^ (near the base of the crust) for B = lO" G. 
Therefore the distribution of currents when the Hall effect be- 
gins to dominate the evolution depends sensitively on the field 
strength. 

4. FIELD EVOLUTION WHEN Qt > 1 

In this section, we first give a qualitative discussion of how 
the Hall effect acts on the magnetic field. We discuss how an 
initial large lengthscale poloidal dipole field generates higher 
wavenumbers, and discuss an interesting time-independent 
solution of the Hall equation. Next, we present short length- 
scale Hall wave solutions for realistic crust models, including 
a first calculation of the elastic response of the crust to the 
wave. We discuss the Hall drift instability found by Rhein- 
hardt & Geppert (2002), and finally discuss the implications 
of our results for GR's non-linear Hall cascade. 

4. 1 . Qualitative discussion of Hall evolution 

In terms of the electron velocity Vg = -J /nee, the induction 
equation ^ is 

— = Vx{veXB). (35) 
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Although this shows a formal relation to the usual MHD ad- 
vection equation, it differs because the electron velocity de- 
pends directly on B. Nonetheless, it is useful to think of 
the field moving with the electron fluid in the same way as 
the field is tied to fluid elements in MHD. It is instructive to 
rewrite equation ( l35t as 



DB „ Ve r 

= B-Vve-B— 

Dt H, 



(36) 



where D/Df = 9/9f + We ■ V is an advective derivative follow- 
ing the electrons, v^, ^ is the radial component of the electron 
velocity, and He is the electron density scale height. The 
first term on the right hand side of equation ( I36> represents 
shearing of the field by the electron flow; the second term 
arises because V- We is non-zero for a vertical current because 
of the background electron density gradient: V- J = gives 
V- We = -Ve ■ '^rie/ He- For example, as currents flow to deeper 
regions, the electron density increases, the electron velocity 
decreases, and the field is amplified by compression. 

As a first example, consider an axisymmetric poloidal field 
Bp = Br(r,9)er + Bg{r,9)ee, supported by toroidal currents 
which we write in terms of the electron angular velocity fte 
as Ve = fleir,9)rsin0e^. In this case, equation i35\ gives 
dBp/dt = 0, and 



dt 



= rsme(Bp-V)ne 



(37) 



showing that any differential rotation of the electron flow 
along the poloidal field lines shears the field, generating a 
toroidal component and a net twist. This is a direct ana- 
log of the winding up of field lines by differential rotation 
(e.g. Mestel 1999). Figure [Tot a) shows an initial dipole field 
and associated toroidal current. In this case, shear in the elec- 
tron velocity at mid-latitudes generates a toroidal field shown 
in Figure [Tof b). The net effect is to twist the original field 
in one direction in the northern hemisphere and the opposite 
direction in the southern hemisphere. The subsequent evo- 
lution is determined by the poloidal currents associated with 
Bcf,. They act back on the original poloidal field, "pinching" it 
at the equator, and generating a quadrupole component of the 
poloidal field, as shown in Figure [TUl c). In this way, the pro- 
cess continues, and the smaller and smaller scale components 
of the field develop. 

As a second example, we consider an axisymmetric toroidal 
field Bt = B^{r,0)e^, supported by poloidal currents Ve = 
Vr(r,9)er + ve(r,6)ee- In this case, equation (I35> gives 



dBj) Ve r I Ve 

-^=-Ve-VB^-B^j^+B^\ 



gCOtt 



(38) 



and dBp /dt = 0, so that an initially toroidal field remains 
toroidal. The terms on the right hand side of equation ( I38> 
describe advection by poloidal currents, divergence of the ve- 
locity flow due to the electron density gradient, and spheri- 
cal g eometry effects. Vainshtein et al. (2000) solved equation 
( I38f for a plane-parallel geometry, in which case it reduces 
to Burger's equation, since only the term involving He sur- 
vives. For a toroidal field with opposite sign s in the northern 
and southern hemispheres, as in Figure flOt b). they showed 
that the toroidal field strength grows at the equator, where the 
electrons flow inwards to higher density, and is reduced in 
strength towards the pole, where the electrons flow outwards 
to lower density (or vice versa depending on the overall field 
direction). They argued that if the field propagates towards 



(b) 




Fig. 10. — The first steps in the evolution of an initial dipole field by 
the Hall effect, (a) The toroidal current supporting the dipole field shears 
the field at mid-latitudes, (b) The resulting toroidal field and associated 
poloidal currents. The net effect is to "twist" the original dipole field, (c) 
The poloidal currents "pinch" the original poloidal dipole field, generating a 
poloidal quadrupole. 



the equator, a thin current sheet will form, rapidly dissipating 
the magnetic energy. It is not clear whether this effect persists 
in spherical geometry and including a poloidal component of 
the field. 

These examples show that the field geometry c an g ive a 
large correction to the Hall timescale estimate of ^2.11 For 
example, in Figure ^| the initial evolution of the field is by 
shearing Z?,. at mid-latitudes. If the magnetic field is con- 
fined to the crust, we expect B,- k, Bg{H /R) to vary across 
a lengthscale H, and 7^ « (c/AnXBg/H). This gives feaii = 
ATTHeeH^ /cBg, so that the horizontal field is th e appropriate 
field to use in the Hall time estimates of ^2.11 If the mag- 
netic field penetrates into the core, then we expect B,- « Bg, 
and Jcf, « (c/4TT)(Bg/R) giving a Hall timescale a factor R/H 
longer This is directly analogous to the sensitivity of the low- 
est order Ohmic decay mode to the boundary condition at the 
base that we found in ^ 

4.2. Time-Independent Solutions 

Not all magnetic field configurations evolve under the Hall 
effect. The existence of time-independent solutions may be 
important for two reasons. First, it implies that dissipation by 
the Hall effect is not necessarily able to cause the field to de- 
crease to arbitrarily small values (assuming that the role of the 
Hall effect is to drive dissipative processes on small scales). 
Second, time-independent solutions may act as a background 
on which Hall waves propagate. The most obvious time- 
independent solution is a force free field with JxB = 0, or 
J = kB, where k is constant along a given field line (e.g. Mes- 
tel 1999). The field profile in this case may or may not reflect 
lengthscales associated with the crust, depending on the initial 
relaxation of the field when the neutron star is born. 

An intriguing time-independent solution is suggested by 
equation (I37> . For an axisymmetric dipole field, the 
toroidal current is 7^ oc sin6' (e.g. Fig. II Oh), implying that 
fte OC Ve/sinO is constant, or the electrons rotate rigidly on 
spherical shells. Equation (I37> then becomes dB^/dt = 
rsinOBridfle/dr). A simple time-independent solution is 
therefore fte =constant, or rigid rotation of the electronsl If 
the electrons rigidly rotate, no shearing of the poloidal field 
occurs, and therefore no evolution due to the Hall effect. This 
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solution implies oc rie, and therefore large currents running 
at the base compared to the top of the crust. If this current 
is sufficiently large compared to currents in the core, then 
deep in the crust the field is Bg ~ -{AttJ ^ ^He / c){n^ I n^.b), 
Br ~ 2Bg(Hg/r), where 7^ and riej, are the current and elec- 
tron density at the base. This solution is not force free, but 
implies an increasing magnetic force with depth; the magnetic 
stresses are poloidal however, and hence may be supported in 
the fluid above and below the crust. 

It is interesting to speculate whether further time- 
independent solutions exist. Equation ( I37> suggests, anal- 
ogous to Ferraro's law of isorotation (Ferraro 1937; Mestel 
1999), further steady-states in which the electrons rigidly ro- 
tate along each poloidal field line. However, the difference 
between the Hall effect and the MHD case is that the electron 
velocity depends directly on the magnetic field. The poloidal 
dipole field is special because the currents are such that no dif- 
ferential rotation of the electrons occurs on spherical shells. 
This is not true in general, and therefore the dipole case may 
be the only steady-state solution of this kind. A toroidal field 
in plane-parallel geometry B = B{x,z)ey (where we assume the 
field is "axisymmetric", d/dy = 0), and with constant n^, does 
not evolve due to the Hall term, since V X ( J X i?) = 0. How- 
ever, Vainshtein et al. (2000) noted that a gradient in n,, gives 
rise to compression of the field (second term on the right hand 
side of eq. 1381 ). and concluded that the field is not steady- 
state in general. This is likely also true once the extra terms 
arising from spherical effects (eq. 1381 ) are included, unless 
the angular distribution of the extra terms can be arranged to 
match (and cancel) the compression term. 

The time independent solution is stable to small perturba- 
tions. In order to see this, let B and J be the magnetic field 
and current density of the time independent state, and 6B and 
6 J be the corresponding perturbations of this state. Lineariz- 
ing eq. Q with respect to the perturbations gives 



ddB 



-Vx 



JxSB + SJxB 



(39) 



Taking the scalar product of eq. ( I39> with 5B /An and inte- 
grating over the volume of the crust gives an equation for the 
magnetic energy in the perturbation 

-dV= I (5JxSB)dV. (40) 



d_ 
dt 



87r 



In deriving eq. ( I40> . we have integrated the right hand side 
once by parts, used Ampere's Law, and assumed the perturba- 
tions vanish on the boundary. Next, we rewrite SJxSB/c as 
the divergence of the Maxwell tensor 



SJx6B 



= V-T, 



where 



T = 



SBSB JB^ 



(41) 



(42) 



A-TT 87r 

We need only the T]^ components of eq. ( I42> because J has 

only a (j> component; J /en,, = -(f)Q,i,rsm6, with 17^ a constant. 
Writing out the integrand on the RHS of eq. ( I40> . using eqs. 
J4U and J42t . and applying Gauss's Theorem, we find that the 
integral vanishes identically. This shows that in the absence of 
any dissipation, small perturbations to the steady state neither 
grow nor decay. Although the realizability of this steady state 
is open to question, the demonstration that it is stable certainly 
bolsters its case for relevance. 



Lastly, we mention that this background solution will decay 
with time at a rate given roughly by the Ohmic time at the base 
of the crust (section l2!2t . 

4.3. Hall waves 

To make further progress, we consider the evolution of 
perturbations to a background magnetic field. GR showed 
that such perturbations propagate as circularly-polarized Hall 
waves. For constant electron density and a uniform magnetic 
field, the induction equation describing plane-wave perturba- 
tions 6B (X exp(/fc • x) is 

dSB ck-B 

— = - (kxSB), (43) 

Ot Anriee 

which may be derived by perturbing equation Q (see also 
Thompson & Duncan 1996). The dispersion relation is 

ck\k-B\ 



47rn-,e 



■ i^k\k-eB\ 



(44) 



where we define u = cB/Airriee = uUciSf oc B/rte- GR discussed 
the suggestion by Jones (1988) that Hall waves could trans- 
port energy to lower densities in the crust, where Ohmic dis- 
sipation is more efficient. They pointed out that as the wave 
moves to lower Wp, the vertical wavelength must increase (for 
fixed uj; see eq. 1441 ). implying a reflection or turning point 
for the wave when kM ~ 1 . In this section, we calculate the 
properties of Hall waves in a realistic crust model. We deter- 
mine the location of the turning point as a function of wave- 
length. In addition, we calculate the strain induced by the Hall 
wave in the crust, and discuss unstable modes in the presence 
of a background electron velocity gradient. 

Force balance in the solid 

The Hall wave in the crust relies on the fact that the solid 
can adjust to accommodate the JxB forces that arise as the 
wave propagates (see also Franco et.al. (2000) for the case 
of magnetic stresses induced by spin down). To see this, 
consider perturbations oc exp(-iwf + /fc • x) of a uniform back- 
ground field with constant ng as before, but now include the 
shear stress (with constant shear modulus /i) and the ion iner- 
tia. Since we assume constant n^, the orientation of the field is 
not important. The induction equation ( I43t has an extra term 
V X {5vX B), giving 

ck ■ B 

-iLuSB=- (kxSB)+kx(uji,xB), (45) 

4-TTnee 

where we write the velocity in terms of the fluid displacement 
as Sv = iuj^. The momentum equation is 



2^ 6JxB 
-puj ^= fik I, 



(46) 



with the terms representing the ion inertia, magnetic force, 
and shear stress. Solving for oj, we find the dispersion relation 



■"Hall ' 



(47) 



where = k^fi/ p and uj^ = {k- B)^ /Airp define the shear and 
Alfven wave frequencies, and tJHaii = ck{k ■ B)/ATmee is the 
Hall wave frequency. 

In the solid crust, which has a finite shear modulus loI > 0, 
there is a low frequency (tj^ <ti u!^,ojs) solution to equation 

ED, 

(48) 



■ 



l+B^/Anp' 
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ForB^/47r <C /i, we recover the usual Hall wave: the induction 
and momentum equations are uncoupled. For B^/4tt ^ fi, 
the emf induced by the fluid displacements is large enough to 
counteract the Hall emf, reducing the Hall wave frequency. 
For typical pulsar fields, this is an important effect only at 
the very top of the crust; for magnetar-strength fields, it is 
an important effect throughout much of the crust. The shear 
modulus is /i « 0.1n,(Ze)^/fl (Strohmayer et al. 1991), where 
a = (3/47rn,)''''' is the interion spacing^, giving B~/4-7rfj, = 
5 X 10"^Bf2(p9ye)"'*/^(Z/26)"2/l In this paper, we assume 
<C 4-Trfi, so that the induction and momentum equations 
are uncoupled in the crust. 

In the fluid above or below the crust, the shear modulus 
is zero, so that the JxB force must be balanced by the 
ion inertia, giving a natural wave frequency u! ^ uja^ ^Haii- 
One might expect that the Hall waves would be almost com- 
pletely reflected at the crust/fluid boundaries because of this 
extremely large mismatch in wave frequency. However, deriv- 
ing the behavior of the strain 9 = d^/dzneai the boundaries re- 
quires a careful treatment of the boundary conditions. In this 
paper, we calculate the response of the crust for the simplest 
case of a vertically propagating wave on a vertical background 
field. We discuss the behavior of the wave at the boundary in 
the Appendix. The key point is that the Hall term enforces 
continuity of the magnetic perturbation across the boundary. 
Therefore, the magnetic stress cannot abruptly change to ac- 
commodate the loss of shear stress on moving into the fluid. 
This implies that ij.9 « BSB/4tt « at the fluid/soUd interface, 
which gives complete reflection. 

Vertically-propagating waves including the density gradient 

We now include the electron density gradient, and study 
vertically-propagating Hall waves in a uniform vertical back- 
ground field. In terms of helicity eigenstates 5B±{z) = 
SBjciz) ± iSByiz), the perturbed induction equation is 



^^^^ 



(49) 



Since V- (55 = 0, 6B~ = constant for these modes. 

We integrate equation j49> through the crust, using a shoot- 
ing method to find the eigenvalues uj for which 5B = Oat each 
fluid/solid interface. Figures II ItHjl show the eigenmodes 
with « = 1, 10, and 100 nodes, normalized so that the maxi- 
mum perturbation is SB = 10'^ G. 

The properties of the eigenmodes can be simply understood 
with a WKB analysis. A good fit to the electron density profile 
in the crust is ne(z) — ne b{z/^Hh)*, with ne t, = 3 x lO-' cm"^ 
and Hi, ~ 250 m the electron density and scale height at the 
base of the crust. Using this to quantize the number of nodes 
n over the crust, J k^dz = mi gives the wave period 



Pn. 



all ■ 



fHaii.b 10^ yrs 



(50) 



,i2 Bxin^ ' 

The WKB solution gives SB oc k^l^ oc nj'^ oc z, in good agree- 
ment with the numerical solutions. The wave reaches a turn- 
ing point when J k^dz ^ 1, or at a depth Ztum ~ ItHb/n^^^. In 
terms of density, the turning point is 



2 X 10'^ g cm" 



n 



4/3 



(51) 




10'» 10>> 

pY^ (g cm-3) 

Fig. 1 1 . — Hall wave with no nodes in a constant vertical background field 
ofB = 10'- G. This wave has a period Pi = 4.4 X 10* yr B^^. We normalize 
the wave so that it has a maximum amplitude SB/B = \. At the top of the 
crust (not shown), 9 drops rapidly to zero (see text). 




10'» 10" 
pY^ (g cm-=) 

Fig. 12. — As Figure fm but for 10 nodes. This wave has a period P\q - 
5.7 X lO" yrBJ^. 



* A convenient way to write this when relativistic degenerate electrons 
provide the pressure is ^/Pe f» 10"^ (Z/26)^'''. 



which requires n k, 300 for ptum ~ lO** g cm""*, or n « 3 x 10^ 
for ptiirn ~ 10'' g cm"^. Above the turning point, SB decreases 
to match the boundary condition SB = 0. Since v oc equa- 
tion ( I49l l gives a steep dropoff SB (xz^, in agreement with our 
numerical results. 
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10'» 10" lO'' 

pY^ (g cm-3) 

Fig. 13. — As Figure fTTI but for 100 nodes. This wave has a period 
/'loo = 580yrB7]. 



The Ohmic dissipation rate for the mode is 



j= dV 




(52) 



where Vg,r is the vertical group velocity, and in the second ex- 
pression we have used constant wave energy flux in the prop- 
agation cavity, and the integral extends between the upper 
turning point Ztp and the base of the crust Zb- This expres- 
sion is simply interpreted as the local damping rate weighted 
by the time spent at that depth. For A:. ^ k±, Vg,, oc fcr' 
so the ratio zk^/p^^^ determines the depth which dominates 

1 /2 

Ohmic decay (for impurity scattering). If k, cx n^,' , the factor 

nl^^ / p^l^ is dominated by the largest densities. We expect that 
up to an order unity prefactor, the Ohmic decay time for the 
mode is set by ilr at the base of the crust, 7"' w ^Haii(f^T)i, 
where (ilr);, w B12/6 when impurity scattering dominates 
(eq. El). 

The elastic response of the crust to the wave 

The strain in the crust may be calculated by balancing the 
elastic and magnetic stresses, since ion inertia plays no role 
in the force balance for the low frequency Hall waves. This 
gives 



dz 47r/x 



-5B. 



(53) 



We show the strain as a function of depth in the lower panels 
of Figures ll 1H13I We write the shear modulus (Strohmayer et 

al. 1991) as /X oc w 10-"' erg cm"-' {z/zb)^^^^ , giving 6 oc 
zT^^/^ below the turning point, and 9 oc zT^I^ above the turning 
point. The maximum strain occurs at or above the turning 



point. Using the WKB scaling, we find that at the turning 
point, 

0.ur„ = 3xlO-^B?2«''/''^, (54) 
where 5Bb is the amplitude of the mode at the base. 

More general cases of Hall waves 

For a horizontal background field, or a non-vertically prop- 
agating wave on a vertical background field, the calculation 
of the elastic response of the crust is more complicated, and 
we leave this for future work. Here, we give the WKB scal- 
ings for waves on a uniform horizontal field B = Be^, and set 
ky = for simplicity, 5B = 6B{z)exp(ikxX-iuJt). The perturba- 
tion equation is 



d^5B, 
~d^ 



+ kl5B. 



■1 



= 0, 



(55) 



giving the WKB dispersion relation w = ±vkx{k^ + k^yl^. 
For k, 3> kx, this gives Pnaii ~ tuaii/ink^Hh). The location 
of the turning point is ztum ~ 3Hb/n^^^, or ptumYe ~ 2 x 
10'^ g cnT^/n'^/^. The turning point for the modes is deeper 
in than the vertically-propagating case; we require n^lQ* for 

Pturn W 10'' g Cm"\ 

Non-zero background current: waves and unstable modes 

Before ending this section, we discuss the case of a hori- 
zontal field in more detail. In several recent papers, Rhein- 
hardt, Geppert, and collaborators have discussed a "Hall drift 
instability" (HDI) driven by curvature in the background field 
(Rheinhai-dt & Geppert 2002, hereafter RG02; Geppert & 
Rheinhai-dt 2002; Rheinhardt, Konenkov, & Geppert 2003; 
Geppert, Rheinhardt, & Gil 2003). In this section, we first 
discuss a simple example of the instability to elucidate the 
basic physics. We then discuss the mode spectrum, and the 
relevance of the instability for small scale Hall waves. 

Following RG02, we study the stability of a plane-parallel 
horizontal field that depends only on height z, B = B{z)ex, with 
corresponding current density J = {c/4Tr)(dB/dz)ey. How- 
ever, for simplicity and to highlight the essential physics, we 
consider the Hall terms only. The perturbations are of the 
form SB(z)exp(ikxX + at) (we set ky = since RG02 find that 
these are the fastest growing modes). Since the background 
field does not evolve due to the Hall term, the perturbed in- 
duction equation is aSB = V X {Sve X i?) + V X (We X SB), or 

(a + Ve-V) SB = (B-V)Sve + (SB-V)ve-(6ve-V) B, (56) 

where = -J/n^e is the electron velocity, and SJ = 
(c /4Tr)(yx SB). 

Writing an equation for SB,, we find the generaUzation of 
equation ( I55t for this case is. 



d^SB, 



■■ K5B, 



, dve/dz 
jykl 



(57) 



{iyk^)\ 

The eigenvalue is purely real^: the solution is either a pure 
growing/decaying mode, or oscillatory. Making a WKB ap- 
proximation, d/dz iky, we find (writing k^ = k^ + k^), 

, dvp 



-(vkk,Y-vki-^ 

' dz 



(58) 



' To show this, multiply equation 1571 by the complex conjugate SB*, and 
subtract the complex conjugate of equation 1571 multiplied by SBz- Integrat- 
ing over volume, and assuming the boundary conditions are such that the 
boundary term vanishes, gives (a^ -a^*) J dV \5B.\' /u^k^ = 0. 
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The first term is the usual Hall waves, a = ±wkk/, the second 
term depends on the electron velocity gradient, and leads to 
instability (cr^ > 0) if the electron velocity gradient is suffi- 
ciently negative, and for sufficiently long wavelengths. 



V dz 



(59) 



A more general integral criterion can be derived by multiply- 
ing equation (|^} by SB- and integrating over the volume. 



dV- 



dV {\5B'X + kl\5BA 



dV\5B-\ 



dvg/dz 



(60) 



showing that dv^jdz must be sufficiently negative somewhere 
in the volume to lead to instability (cr^ > 0). We have numer- 
ically integrated equation ( I57> for the simple case of constant 
He and a quadratic dependence of B on depth, confirming the 
scalings of equations ( I58> and ( I59> . 

The basic physics of the instability can be understood by 
stepping through the evolution of an initial vertical perturba- 
tion ^B,, and recalling that the field is advected by the electron 
flow. First consider short wavelengths where the dv^/dz term 
is unimportant. The current perturbation 5Jy = -ik^-6B,c/4TT 
corresponds to a sheared electron flow, which gives rise to a 
perpendicular field component aSBy = -ik^-BSv^.x through the 
first term on the right hand side of equation (I56> . Associ- 
ated with SBy is a current 5J~ = ikySByc/AiT, which in turn 
acts back on 5B^, a6B, = -ik^-BSv^ -, thereby driving the Hall 
wave. When the dvg/dz term dominates, the first stage of 
this evolution is different. Shearing of the initial perturba- 
tion by the background electron flow generates SB^ in the first 
step, aSBy = -SBr{dvg/dz), transferring energy from the back- 
ground field into the perturbation. Equation (I58> shows the 
competition between the internal and external shearing, mea- 
sured by the ratio of Hall wave frequency to the shear fre- 
quency dve/dz- 

The criterion for instability we derive above generalizes that 
of RG02 to include the density variation with depth (see also 
Rheinhardt et al. 2003), and clarifies the nature of the in- 
stability: a background shear in the electron velocity drives 
growth of long wavelength perturbations whose characteristic 
Hall frequencies are less than dv^jdz- Short wavelength Hall 
waves, with frequencies greater than dvg/dz, are unaffected. 
This differs from the conclusions of RG02, who argue that 
HDl may be responsible for the growth of small-scale pertur- 
bations reported in some numerical simulations (e.g. Shaly- 
bkov & Urpin 1997; Hollerbach & Rudiger 2002). Perhaps 
numerical instabilities are rather responsible for the excess 
small scale power*. 

Is the HDI relevant for the evolution of field in the crust? If 
the picture of a turbulent Hall cascade is correct, the instabil- 
ity probably does not change the long-term evolution of the 

* We found such an instability in trying to develop a 3D code to simulate 
the Hall effect. It is straightforward to see this in ID in Cartesian coordinates. 
Using second-order centered differences in space and explicit time-stepping 
(first order, or, using Runga-Kutta, second order accurate), one can analyti- 
cally show that all wavelengths are unstable! The growth rate of the instabil- 
ity for wavenumber kona grid of spacing Aa: scales as oc (kAjcf" , implying 
large lengthscale waves grow only after many box crossings, but grid-scale 
waves have a growth rate of order the time step. Including Ohmic diffusion 
only cures this instability if a prohibitively short time step is used. In ID, the 
instability is easily fixed by using implicit time- stepping, but this cannot be 
simply extended to 3D using operator splitting. 



field, since intermediate scales will "fill in" as the cascade de- 
velops. If the background field varies on a lengthscale H, the 
most unstable modes are those with the largest ky compatible 
with equation ( I59l l. implying ky ^ 1 ///, or corresponding to 
spherical harmonic / of tens. Rheinhardt et al. (2003) then ar- 
gue that an initially dipolar field will develop small latitudinal 
scale structure; however, this conclusion is sensitive to the ini- 
tial distribution of currents in the crust. An additional point is 
that changes in the crustal field may not propagate to the sur- 
face because of the very efficient reflection at the solid/fluid 
interface (see discussion in the Appendix). Therefore predic- 
tions of complex surface fields produced by HDI (Geppert et 
al. 2003) may be premature. 

4.4. Nonlinear Coupling of Oscillation Modes 

We now discuss the implications of our results for the turbu- 
lent Hall cascade suggested by GR. They estimated the non- 
linear transfer time between different modes to be longer than 
the mode period, and proposed that the cascade would be a 
weak cascade. In this case, the expected scaling of magnetic 
energy is B^(k) cx k~^, with a cutoff due to Ohmic dissipation 
at a scale A ~ L/VIt, where L is the scale on which the turbu- 
lence is being driven. However, the simulations of Biskamp 
et al. (1999) show a different scaling, B^{k) oc k'''!^ appro- 
priate for a strong cascade (to derive this scaling, follow the 
argument of GR §4, but choosing their eq. [47] rather than 
eq. [48]). In this case, the Ohmic cutoff is for A L/iyLrfl^. 
The difference arises most likely because of GR's assumption 
of isotropy: a strong cascade implies mode periods longer 
than the non-linear transfer time, which suggests a cascade 
perpendicular to the magnetic field direction, for which the 
Hall wave periods are long. Unfortunately, Biskamp and 
coauthors did not present an analysis of anisotropy in their 
numerical results. Such a study would be extremely interest- 
ing. 

How many modes are involved in the cascade? If the energy 
transfer is limited by Ohmic dissipation, then for a strong cas- 
cade we find 



{nrfl^ « 3 X 10"* 



3/2 



(61) 



where we use the result for impurity scattering for VIt. 
The turning point for the mode with n = n^ax (assum- 
ing vertical propagation and a vertical field) is at p^ax ~ 
10'' g cm-3 {Q-^lBnf. 

In ^4.31 we calculated the maximum strain 0n,ax for a Hall 
wave with given n and amplitude at the base SBj,. As n in- 
creases, ^max increases, implying that the crust will yield for 
large enough n. However, this increase with n is offset by 
the decreasing amplitude 6Bi, with n if there is a turbulent 
cascade: for a strong (weak) cascade, 5B oc n"^/* (cx n"'). 
The value of strain at which the crust will yield 6'yieid is un- 
certain (e.g. Ruderman 1976), but probably lies in the range 
6'yieid ^ 10"^- 10"^. For the weak cascade, 6'tnax decreases more 
slowly with n. Using equation ( I54t for vertically-propagating 
waves, we find for this case that the crust will yield for modes 
with 



n > 6 X 10^ 



^yield ^ ^ j^-ip f SBh.o 
\ B 



To=4 j 



-13 



(62) 



where 5Bi,,o is the amplitude of the mode at the largest length- 
scale, and (3 = 9/4. For a strong cascade, the prefactor is 
2 X 10'' (corresponding to fir w 10*) and (3 = 3.2. Given our 
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estimate for n^ax (eq. 16 II ). we therefore do not expect the 
cascade to be limited by crust yielding, but rather Ohmic dis- 
sipation. However, we stress again that we have not yet cal- 
culated the crust response for arbitrary wave propagation and 
background field directions. Therefore, equation ( I62t should 
be taken as an initial estimate only. 

5. SUMMARY AND CONCLUSIONS 

In this paper, we have addressed various aspects of mag- 
netic field evolution in neutron star crusts. Our aim has been 
to place previous studies of the physics of the Hall effect in 
context: first, by addressing when and for which neutron stars 
the Hall effect is important relative to Ohmic decay; second, 
by calculating the allowed "initial conditions" by following 
the Ohmic decay prior to dominance of the Hall effect; third, 
by discussing for the first time the properties of Hall waves 
in a realistic crust model, and the implications for the Hall 
cascade. 

We now summarise our main results and conclusions: 

(i) The relative importance of Ohmic decay and the Hall 
effect is summarized in Figure |4] which shows the fir = 1 
contour in the B-T plane. In isolated neutron stars with 
B < 10^^ G, the dominance of the Hall effect depends sen- 
sitively on the impurity parameter Q: recent calculations by 
Jones (2001) indicate 2 > 1 for densities greater than neutron 
drip (see also de Blasio 2000 for densities lower than neutron 
drip), in which case the Hall effect never dominates the evo- 
lution of crustal currents, considerably simplifying the evolu- 
tion. If 10"^ however, as in the original Flowers & Ruder- 
man (1977) estimate, Ohmic decay dominates the Hall effect 
only for a time fswitch ~ 2 x 10"* years/Bjj, after which the con- 
ductivity is set by impurity scattering, giving fir « Bn/Q- In 
accreting neutron stars with B < 10'^ G, Ohmic decay always 
dominates, since the crust is heated (T > 10** K, e.g. Brown 
2000) and may also be impure Q > 1 (Schatz et al. 1999). 

(ii) Figure |9l summarizes the evolution of currents in the 
crust while Ohmic decay dominates (t < fswitch)- The electri- 
cal conductivity at these early times is set by phonon scat- 
tering. This figure relates previous work on either Ohmic 
decay eigenmodes, or on self-similar solutions, in a simple 
way. For a particular time t, Figure |9] gives the density 
down to which the currents have significantly decayed. For 
p > pc, no significant evolution of the currents has occurred. 
For p < Pc, the current distribution is self-similar, being dom- 
inated by the Ohmic decay mode with turning point at p= pc- 
Because fswitch l/B^, the allowed distribution of currents 
when the Hall effect begins to dominate is very sensitive to 
B. For B K, 10'"^ G, pc lies at or above the top of the crust; 
for B !v lO" G, Pc lies at the base of the crust. At late times, 
when Pc is greater than the initial location of the currents, the 
current distribution evolves in a self-similar manner and the 
surface field decays as a power law (eg. 1311 ). 

(iii) The basic physics of the Hall effect is that the field in 
the crust is advected by the electron flow. This is summarized 
in Figure^ which considers a variant of the usual laboratory 
demonstration of the Hall effect. Figure^] shows the evolu- 
tion of an initial poloidal dipole field, in which the toroidal 
currents initially "twist" the field; the resulting poloidal cur- 
rents then generate a quadrupole poloidal field via a "pinch". 
These considerations led us to a simple time-independent so- 
lution: a poloidal dipole field, in which the electrons rotate 
rigidly, does not evolve due to the Hall effect. 



(iv) We discussed the properties of Hall waves in the crust. 
The simplest example is a background magnetic field which is 
uniform and vertical (for example, emerging from the core). 
The vertical wavelength of the wave grows towards the top 
of the crust (as pointed out by GR), and we showed that the 
turning point location is ptum ~ 10'^ g cm"-'/n^''-'. Below the 
turning point, the wave is able to propagate and its proper- 
ties are well-described by WKB scalings. Above the turning 
point, the wave evanesces towards the top of the crust. 

(v) A horizontal field (e.g. as may be produced if the cur- 
rents supporting the field are predominantly in the crust) has 
similar wave-like solutions at short wavelengths. However, if 
there is a background current (the background B varies with 
height), long wavelength perturbations become unstable. The 
physics of this instability is that the shear in the background 
electron velocity drives unstable growth, overcoming the in- 
ternal shear of the wave which would otherwise drive an os- 
cillation. Therefore we agree with the recent work of Geppert 
and collaborators that there is an instability, and we have clar- 
ified its basic physics. However, we argued that the relevance 
of this instability for the field evolution in the crust, and par- 
ticularly for the observable surface field, is not clear. 

(vi) We have made a first attempt to address the question 
of whether Hall waves strain the crust beyond its yield point. 
We find that the maximum strain is at or above the turning 
point for the mode (eq. 1541 ). However, our calculation is 
limited to the special case of vertically-propagating waves on 
a vertical background field; we leave further calculations of 
the response of the crust to non-vertical background fields 
and wavevectors for future work. In particular, the force bal- 
ance is important to consider for magnetar strength fields, for 
which B^ > Airp, directly coupling the induction and momen- 
tum equations in the crust (see the recent discussion in Arras 
et al. 2004). 

(vii) Finally, we briefly discussed the non-linear evolution. 
Recent numerical work on Whistler turbulence (Biskamp et 
al. 1999) has demonstrated GR's conjecture that the Hall ef- 
fect would create a "cascade" of energy to small scales. How- 
ever, the scalings of the cascade imply strong, anisotropic tur- 
bulence rather than weak, isotropic turbulence as assumed by 
GR. The Hall wave eigenfunctions calculated in §4 (again, for 
vertically propagating waves and a vertical field) show that 
crust breaking is probably not important for standard radio 
pulsar strength fields. 
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APPENDIX 

BOUNDARY CONDITIONS FOR HALL WAVES 

In this Appendix we discuss the boundary conditions for Hall waves in the crust. First we argue that the disparity in propagation 
speeds between solid and liquid implies near total reflection of wave energy at the boundaries. We present the results of simple 
toy models to quantify the transmission coefficient. We then derive the boundary conditions. Lastly we specialize the boundary 
conditions for vacuum solutions in the liquid. 

For elastic response of the crust, i.e. neglecting irreversible processes such as plastic flow, we naively expect Hall wave energy 
to be almost perfectly trapped in the crust. For a wavelength of order the size of the crust, the Hall wave propagation speed is 
vastly smaller than the Alfven speed, implying a large impedence mismatch between the solid crust, liquid core and liquid ocean. 

To get an analytic handle on the size of the reflection and transmission coefficients, consider a region with uniform density, 
constant vertical background field, crust with constant shear modulus for z < and ocean with zero shear modulus for z > 
(see Blaes et al. 1989 for the Alfven wave case). This toy problem is of limited validity, since the density is not constant, and 
also since in reality the WKB assumption is violated (vertical wavelength in the ocean ^ scale height), but we can gain useful 
intuition without much effort. We consider a vertically propagating Hall wave incident from below, and then include all possible 
reflected and transmitted waves. The two circular polarization states are decoupled in the boundary conditions, hence we need 
only include reflected Hall and elastic-Alfven waves and a transmitted Alfven wave. We use the boundary conditions derived 
below to determine the amplitudes of reflected and transmitted waves. We leave this straightforward calculation as an exercise 
for the reader, and present the results. In the limit <C 47r/x near the top of the crust, we find the amplitudes of the reflected 
Hall, reflected elastic-Alfven, and transmitted Alfven waves are 

5B{^ei Hall) — - f 1 + — — — ) ^fi(inc Hall) 
V I'Alfven / 

AR ~ VHallVAlfven ^p 

CD(ref Alfven) — J ''"(inc Hall) 

^shear 

an Alfven) — 

—5B 

(inc Hall), (Al) 

^Alfven 

yielding an energy transmission coefficient 

T = 2vHaii/vAifven ^2.4 X IQ-'^ [b ,^nP{yv)]~' . (A2) 

Hence Hall waves experience negligible damping due to wave leakage in this simple case. Interestingly, we find that both the 
reflected and transmitted Alfven waves have a very small amplitude compared to the Hall wave. In this simple toy problem, the 
(observable) field in the liquid is composed of the constant background field (generated in the core) as well as the perturbation. 
Since the amplitude of the perturbation is infinitesimally small in the liquid as compared to the solid, the observable external field 
may be somewhat different than the field in the crust. If this behavior persists for more realistic field distributions, the observable 
external field may not be the same as the field in the crust; the neutron star is in effect "hiding" its field. This effect is of course 
larger when the currents are mainly in the crust rather than the core. 

If we relax the restriction of vertical propagation (i.e. allow finite k±), there is a qualitative change in the nature of the waves. 
The low Hall wave frequency implies the reflected elastic-Alfven wave becomes evanescent with e-folding distance ^ k']^, and 
the transmitted Alfven wave has wavevector nearly perpendicular to the background magnetic field. Both polarization states must 
now be included for the reflected elastic-Alfven wave and transmitted Alfven wave. In addition, the two polarizations of reflected 
Hall wave correspond to propagating and evanescent waves. With so many waves included, the analysis is complicated, and a 
nontrivial check of both our numerical results and the boundary conditions derived below was that energy was conserved over 
a large range in frequency and k±. For incident Hall wave of amplitude unity, the numerical results can be summarized as : (i) 
the amplitude of the reflected propagating Hall wave is 1 + O (k± / A;z,Haii) , (ii) the amplitude of all other reflected waves scales as 
O {k±/k2„Haii), (iii) the transmission coeffient is independent of k±, and agrees in order of magnitude with eq lA2l Hence even in 
this more general case the crust still acts as a nearly perfect resonant cavity, and the wave is essentially unobservable since it's 
amplitude in the liquid is so small. 

We now discuss the boundary conditions. We do not linearize the magnetic field, but for simplicity we use the linearized shear 
stress. We use the notation [/] = (/)soiid-(/)iiquid, and the normal to the boundary is in the z direction. By integrating over a small 
region near the solid-liquid boundary we find 

[B,] = from V-B = (A3) 

[J,]=0 fromV-J = (A4) 

[^,] = from finiteness of V- 1 (A5) 

(/:- 2^/3) V- 14 + /^ + 4,,) + (47r)"' (^BiB,-^B^Si?j +g^z[p]Siz = from momentum equation (A6) 

= from induction equation. (A7) 

The last term in equation iAbi ar ises from the buoyancy force -e^gSp = e^g'V- (p^;), and we have ignored the perturbed gravi- 
tational potential. Equation iA5i is required since V-^ would be infinite if was discontinuous across an infinitesimally thin 
layer 
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The Hall effect, like Ohmic diffusion, involves a higher spatial derivative than the ion advection term in the inducti on e quation. 
Hence an additional boundary condition comes from the requirement that these terms remain finite in equation (IA7> . If the 
magnetic field was discontinuous across the effectively infinitely thin boundary, the Hall and Ohmic electric fields would diverge. 
Hence we must add the boundary condition 

[B±] = from finite Hall voltage. (A8) 

This boundary condition is similar to the no-slip condition from viscous hydrodynamics^. 

The hor izon tal stress boundary conditions can be greatly simplified to give a useful result for the strain. Taken together, 
equations (IA3> and (IA8> imply [B] = 0, which has the interesting consequence that the magnetic stress drops out of the boundary 
condition in ea lA6l Since ^ = in the liquid, the horizontal components of the (symmetrized) strain must vanish on the solid side 
of the boundary, 

(6,z + ev),oiid = for i = x,y. (A9) 

As the magnetic field in the crust evolves by the Hall effect, equation ( IA9I I is not satisfied in general over the entire crust. As 
the wave approaches the boundary, it must adjust in order to satisfy this condition. On qualitative grounds, we expect this would 
reflect the magnetic energy even if the wave is not short lengthscale. This result suggests that even for the fully nonlinear case, // 
the currents are generated in the crust there may be a significant difference between crustal field and observable external field. 

As discussed by Kingsep et al. (1990), the higher spatial derivatives give rise to a boundary layer effect in Hall-MHD (without 
the elastic shear stress), with ion advection dominating on large scales and the Hall effect on small scales (shorter than the ion 
penetration depth). As a test of the boundary conditions, we have redone the transmission-reflection problem of Blaes et al. (1989) 
for Alfen waves to escape from the crust, now including the Hall term in the induction equation. We find that for typical Alfven 
wave frequencies we recover the Blaes et al. (1989) Alfven-wave transmission coefficient, consistent with the idea that the Hall 
term only affects Alfven waves quite close to the boundary, but cannot alter a global property such as transmission out of the 
crust. 

Lastly, we note that if the currents are small in the liquid (say due to low conductivity), or equivalently if the field is near the 
vacuum solution, one finds an interesting boundary condition on the field in the solid. Since 7, = in the liquid, and [J^] = 0, 
7; cx kxBy-kyBx = In the solid also. We can always choose = without loss of generality, again implying By = in the solid 
at the boundary. 
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